Anomalous diffusion for a correlated process with long jumps 



Tomasz Srokowski 

Institute of Nuclear Physics, Polish Academy of Sciences, PL - 31-342 Krakow, Poland 

(Dated: June 21, 2011) 

We discuss diffusion properties of a dynamical system, which is characterised by long-tail distri- 
butions and finite correlations. The particle velocity has the stable Levy distribution; it is assumed 
as a jumping process (the kangaroo process) with a variable jumping rate. Both the exponential 
and the algebraic form of the covariance - defined for the truncated distribution - are considered. 
It is demonstrated by numerical calculations that the stationary solution of the master equation for 
the case of power-law correlations decays with time, but a simple modification of the process makes 
CN| , the tails stable. The main result of the paper is a finding that - in contrast to the velocity fluctu- 

ations - the position variance may be finite. It rises with time faster than linearly: the diffusion 
is anomalously enhanced. On the other hand, a process which follows from a superposition of the 
_ Ornstein-Uhlenbeck-Levy processes always leads to position distributions with a divergent variance 

. which means accelerated diffusion. 

(N 

Keywords: Diffusion; Jumping process; Correlations; Stable Levy distribution 

o 

dJ ; I. INTRODUCTION 

Stochastic force in the dynamical description of a physical system is usually assumed as a white, Gaussian noise. 
C$ ■ That approximation may be sufficient if the correlation time is short, compared to the time-scale of the system. 
However, correlations of the stochastic force are important in many physical problems. Slowly decaying correlations 
. are present in various phenomena including chemical reactions in solutions [H, ligand migration in biomolecules 0, 
atomic diffusion through a periodic lattice Q , spectroscopy and hydrodynamics [H, Long-time memory may 
also occur in climate, physiology, financial markets and earthquake research The paper Ref.[l(| reported that 

the normalised prices of both Babylonian and English agricultural commodities exhibited persistent correlations of 
a power-law type over long periods of up to several centuries. Moreover, long-time correlations emerge in effective, 
low- dimensional stochastic descriptions of complicated deterministic systems; they are produced by the procedure of 
q \ fast degrees of freedom removal [ll], [ijj . Long correlations of the velocity are often connected with the anomalous 
1-1 . diffusion: the position variance rises either slower or faster than linearly with time. Examples of anomalous diffusion 
are frequently encountered in disordered media |13j . It has been demonstrated recently that the subdiffusive motion 
of a polymer is related to scaling of the monomer mean square displacement and the velocity autocorrelation function 
' is algebraic [l4j • The importance of strong memory (long autocorrelations) - as well as long tails of the distribution- 
Q\ , in fluctuations of company profit was stressed in Ref. |15| : simple, uncorrelated models are a poor approximation 
to real market phenomena. Systems with finite correlation time of the stochastic force may be described by the 
integro-differential Langevin equation. If the memory kernel has the same form as the force autocorrelation function, 
\Q • fluctuation-dissipation relations are satisfied (16|. 

Anomalous transport can be described by jumping processes. Continuous time random walk theory (CTRW) 17], 
in particular, predicts subdiffusion which results from the non-Poissonian form of the waiting-time distribution and 
makes the process non-Mar kovian. That property is mathematically expressed as a fractional time-derivative in the 
, Fokker-Planck equation. Anomalous diffusion, both super- and subdiffusion, is predicted by a Markovian version 
of CTRW if one allows for a variable jumping rate (is) . Besides the waiting-time distribution, CTRW involves the 
jumping-sizc distribution. Due to the central limit theorem, stable distributions are distinguished. They can have 
either the form of a Gaussian or a general stable Levy distribution, which is characterised by long tails ~ |x| _1_a , 
where < a < 2 is the stability index (Levy nights) [TiJ |2(j. Then, the Fokker-Planck equation contains the 
fractional derivative in respect to the process variable. The above asymptotics implies that variance, as well as all 
higher moments, must be divergent for any time (the accelerated diffusion). This property means that infinite jumps 
are performed during a finite time; for that reason Levy flights are unphysical in dynamical problems. The diffic ulty 
can be avoided if one introduces the Levy walk for which the jump size is strictly related to the time interval [17j . 
Those systems can exhibit anomalous diffusion. This is the case for the transport of two-level atoms in optical molasses 
derived from counterpropagating laser beams [2lj | . The atom which stays for a long time in the potential well may 
accumulate enough energy, due to spontaneous emissions, to jump above the barrier and perform a long flight. If the 
depth of the optical potential is small enough, both the momentum autocorrelation function and the position variance 
obey the power-law dependence; diffusion may be anomalous. One can expect that for the Gaussian processes the 
system is close to the thermal equilibrium; the detailed balance is then satisfied. It has been recently demonstrated 
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that CTRW with a constant driving force satisfies the Einstein relation for an arbitrary memory kernel [22} . 

A jumping process, which is well suited to modelling long-time correlations, is called the 'kangaroo process' (KP) p3j . 
It is defined by the Poissonian waiting-time distribution with a variable rate. Instead of the jump-size distribution, KP 
assumes the distribution of the process value just after the jump. As a result, it cannot be described by a differential 
- either second order or fractional - Fokker-Planck equation, and then it is well suited as a model of non-diffusive 
processes with strong collisions. KP is a natural formalism to handle the fully developed boundary-layer turbulence 
since a local description is inadequate in this case [24j ; the velocity profile for the Newtonian shear flow is governed by 
the integro-differential equation of KP. Moreover, a stochastic description of the hyperfine structure of spectral lines 
is possible by KP since correlations of the radiation emission/absorption operator can be easily taken into account 
[25l | . For the same reason, KP was applied to describe the stochastic Stark broadening, which involves correlations 
~ 1/t 4] . It can also serve as a model of coloured noises in dynamical problems [26l — f28l |\ 

Considerations involving the correlation problem require some care for the case a < 2 since the covariance function 
is infinite. One can introduce a simple modification of the process by cutting the distribution at some large value. 
Such a truncation procedure is realistic since all physical systems are finite. For example, truncation of the velocity 
distribution emerges in the natural way in optical molasses since the Doppler force becomes large at high momentum 
[2lj . Processes with truncated distributions have finite (co— )variance, besides that their properties do not differ from 
those for the stable distributions [29|, |3(| • 

The aim of this paper is to demonstrate that diffusion - in the presence of the Levy flights - qualitatively depends 
on both the correlation form and specific choice of the noise. We consider two models of the particle velocity: a 
jumping process (KP) and the Ornstein-Uhlenbeck-Levy (OUL) process. The paper is organized as follows. In Sec. II 
KP is defined; its properties and stationary solutions are discussed. The distributions of the position and diffusion 
properties of the system are discussed in Sec. III. In Sec. IV the results are compared with the predictions of a model 
which is based on the OUL process. 



II. THE JUMPING PROCESS 



The process is defined in terms of two probability distributions Q(£) and Pp(r). Q(£) determines the step-wise 
process value just after a jump. More precisely, the probability that one jump occurred in an infinitesimal time 
interval (0, At) is given by v(^)At, where i^(£) denotes a given, variable jumping rate. Therefore, the waiting-time 
density distribution, Pp(r), is Poissonian, 

Pp(t) = K£)e-^ )r ; (1) 

it determines the probability that the jump occurs for the first time in the interval (t, t + At). Immediately after the 
jump, the probability density of £ becomes Q(£) and the process value remains constant until the next jump. The 
process is stationary and Markovian, and it can be described by the infinitesimal transition probability £' — > £: 

p tr (£, At^, 0) = [1 - i/(£')At]5(£' - f ) + v(OAtQ(0- (2) 

Since the time interval At is small, only one jump is taken into account in the above expression. The first term in 
Eq.@ corresponds to the event that no jump occurred. From Eq.©, the master equation directly follows: 

-p(£, t) = -i/(0p(£, t) + Q(0 j (3) 
The stationary solution of Eq. ^ is given by 

where the frequency z/(£) is averaged over Pi(£). Existence of the integral in Eq.Q imposes appropriate requirements 
on the functions Q(£) and However, Pi(£) is not the only stationary solution of Eq.©: if limj_^ a f(£)/(£ — a) 

is finite for a constant a, a delta function also constitutes a solution. Then, for a = 0, Eq.Q implies 

m) = m (5) 



Therefore, the master equation © can have the regular solution (Eq.Q), the singular one (Eq.©) or both, according 
to the specific choice of the functions Q(£) and 
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The process is characterised by a broad spectrum of possible covariance functions. We define that quantity in the 
stationary limit by a joint probability P as 

c(t) = ma(t+T))/m 2 ) = J J^p(Lt+r;e,t)d^e/m 2 ) = J J^'piu+rie^^dc/m 2 ), (6) 

where we assume that all the above quantities exist and are finite. C(t) can be evaluated by summing up over all 
possible sequences of jumps [Hj]. If both P\(£,) and are even functions, Eq.© takes a simple form 

/>oo 

C(t) = 2 £ 2 exp(-^)t)Pi(£)^/(? 2 )- (7) 
Jo 

The integral in Eq.0 can be expressed in the form f^? \ £ 2 exp(— t / (^)t)Pi(^)(d^/dv)du which is essentially a Laplace 

transform. We want to construct the jumping process for a given covariance by an appropriate choice of the (even) 
function This can be achieved by inverting the Laplace transform from C(t); then follows from a simple 

differential equation 

§=2e 2 Pi(0/CM, (8) 

where C{v) denotes the inverse Laplace transform of C{t). u(x), calculated from Eq.®, defines the required process 
as soon as Q(x) is an even function. In particular, for the power-law covariance 

C(t) ~ 29r(6)t- e {6 > 0), (9) 

Eq.© yields 
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(10) 



In this paper, we assume the stationary distribution Pi(£) in the form of the stable, symmetric Levy distribution. 
It is given by the following Fourier transform 

1 f°° 

Pi(0 = ~ cxp{-Dk a )cos(k£)dk, (11) 
n Jo 

where the order parameter a G (0, 2]. Except for the case a — 2, which corresponds to the normal distribution, Pi(£) 
possesses the algebraic tail ~ Then the variance is divergent and the covariance definition, Eq.®, does 

not directly apply. To overcome that difficulty, one introduces modifications to that definition, e.g. the codifference 
[32l l33l] . which can serve as an estimation of the memory loss, or introduces an infinite cascade of Poissonian correlation 
functions [34| . On the other hand, we can approximate the Levy distribution by introducing a cut-off at some large 
value of the argument. Those truncated distributions |3 |. [3 5j-l39f are realistic in respect to physical applications, and 
they agree with the exact ones up to arbitrary large £ [29j; convergence to the normal distribution is very slow. In 
the following, C(t) will be understood in the sense of the truncated distributions. 

In the simplest case of constant jumping rate, = Vq =const, the master equation, Eq.([3]), takes the form 

~§iP(£i t) = v o[Q(€) ~ *)]• Its solution with the initial condition 0) = <5(£), 

p(C, t) = 5{^- Vot + 0(0(1 - *~ Vot )' (I 2 ) 

converges with time to the unique stationary distribution Pi(£) = Q{Ci- The exponential form of the covariance 
follows directly from Eq. ([7]). 

Equality of the distributions Q(£) and Pi(0 does not hold for the variable jumping rate cf. Eq.Q. We 

consider the case of the power-law covariance, Eq.Q. To ensure the existence of Pl(£)j we impose a condition on 
the parameters a{6 + 1) > 2; then (v(£)) < 00. On the other hand, the singular solution Pa(£) (Eq.©) also exists if 
9 < 3 because v(£) ~ \£\ 3 ^ 8 for |^| -c 1. Which of those two solutions, Pi(£) or Pa(^)) i s stable? This problem was 
considered in Ref.[4(| for the constant Pi(£) on a finite interval. If one approximates the interval length r in Eq.([T]) 
by its average, the process is fully determined by the distribution of r (P(t)) which, in turn, uniquely follows 

from Pl(£)- However, we are interested in the process value at a given time t, i.e. the last interval in the evolution 
must contain t. That requirement imposes a bias, namely longer intervals are more probable. As a result, the effective 
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FIG. 1: (Colour online) Time evolution of the density distribution p((.,t) for algebraic C(t). Solid lines correspond to the 
following times: 0.05, 0.5, 2, fO and 40 (from top to bottom at large £). The distribution Pi(£) is marked by the red dashed line 
and Q(£) by green dots. Inset: relative decay rate of p(£,t) (dots) and estimation of that rate by the function f(t), according 
to Eq.([I5| (solid line). 



P(r) is modified: it acquires a time dependence and the effective variance dwindles with time to zero. Therefore, 
the solution Pi(£) appears unstable, and it decays to the delta function. In the following, we will demonstrate by a 
numerical analysis that Pl(£) in the Levy form is also unstable. 

Simulation of the stochastic trajectory requires sampling of the process value £ from the distribution Q(£), which 
is given by Eq.Q. Details of that procedure - for the Levy distributed Pi(£) and the power-law covariance - are 
presented in Appendix. The procedure is relatively simple if 9 = 1 and a — 3/2. In the following, we restrict the 
numerical simulations to that case. Consecutive time intervals are given by the distribution JTJ). The time evolution 
of the density distribution p(£,t) is presented in Fig.l, where the initial condition p(£, 0) = Q(£) was assumed. The 
tail quickly converges to the Levy shape |^| _1 ~ a but its relative strength falls with time. The rate of the decay can 
be estimated by a function which consists of two power-law segments; it is given by 



/(*) 



d*- ft for t < 7.5 
c 2 t-^ for t > 7.5, 



(13) 



where /3j = 0.584, 02 — 1.03, c\ — 0.586 and C2 = 1.46. Decay of the tail is compensated by an increase of the 
probability density near £ = 0: we observe convergence to Pa(£)' 

The nonsingular solution becomes stable for large |£| after an appropriate variable transformation, £' = <?(t)£. The 
function g(t) follows from the requirement that the tail of the new distribution, P{(£') ~ a_1 i is time-independent. 
Using the identity P{(£') = Pi(£)RM'|, where P 1 (£) - /(*)|^"| _at_1 , yields 



9(t) 



c'^'i for t < 7.5 
d,t& for t > 7.5, 



(14) 



where /3- = Pi /a and c- — 



-1/q 



Numerical calculations, presented in Fig. 2, confirm that the tail of the limit 



distribution is stable and coincides with the Levy one. The decay of the distribution Pi (£) with time accelerates the 



decline of the covariance (El): C 



(i) 



e// (r) = /(r)C(r). On the other hand, C^(r) = (£'(t)£'(t H 
Both covariance functions are presented in Fig. 2; they follow the dependence r _1_/31 and r 



9(t)C(t). 
respectively. 



III. DIFFUSION 



Solution of the Fokker-Planck equation for the standard theory of the Brownian motion predicts that the position 
variance rises linearly with time. Deviations from that pattern may be caused by nonhomogeneity of the medium 
[l3l [l8[ , the presence of regular structures in the phase space (4lT | or memory in the system (l~7| . In general, we observe 
the subdiffusion and superdiffusion (the enhanced diffusion), when the variance rises slower and faster than linearly 
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FIG. 2: (Colour online) Distribution of KP values as a function of the transformed variable f = g(t)£. Inset: effective 
autocorrelation functions C^L(t), calculated from the definition at t = 1. The distribution Q(£) was truncated at £ = 10 4 . 
The solid lines have the slopes 1.584 and 1.195. 

with time, respectively. Moreover, the accelerated diffusion corresponds to the case of infinite variance for any time. 
If the process is stationary and the velocity autocorrelation function C(t) exists, the time-dependence of the variance 
is directly related to C(t) by a simple identity 



where the average is performed over the stationary distribution. The above formula means that slow decrease of 
C(t) implies fast diffusion. The non-Markovian CTRW with the Gaussian jump-size distribution predicts subdiffusion 
which results from long waiting times inside the traps. Anomalous diffusion may originate not only from temporal 
characteristics of the system but also from the spatial structure of the medium. It is the case for transport on fractals 
[42| and for the Markovian CTRW with position-dependent jumping rate (43|; all kinds of diffusion are present in 
those systems [T^]. The stochastic, additive driving force in the Langevin equation with the stable and non-Gaussian 
distribution leads to the accelerated diffusion. The same conclusion holds for CTRW with the stable jump-size 
distribution but introducing a power- law truncation may result in heavy tails with convergent variance (38j . 
We look for the distribution p x (x,t), where the variable x is given by the equation 



and v(t) follows from the jumping process. In the case z/(£) =const, which corresponds to the exponential covariance, 
we assume v(t) — The distributions p x (x, t) were calculated from individual trajectories x(t); consecutive process 
values were sampled from the Levy distribution by a standard procedure |44|. Fig. 3 presents the results as a function 
of the Levy parameter a. The tails of all distributions for a < 2 have the shape which implies infinite variance 

and accelerated diffusion. The apparent width of p x (x,t) diminishes with a and rises with time. 

The case of variable jumping rate involves a coupling between the process values £ and the waiting times r, cf. 
Eq. (|10p . Since a large £ corresponds to a small time interval r, its influence on the particle displacement is limited, 
and the resulting variance may actually be finite. The condition for the convergence of the variance can be estimated 
in the following way. We assume v(t) = £(t) and take into account only the tail of the Levy distribution, ~ ^ _Q_1 . 
Then Eq. lfTt))) yields z/(£) <~ £( 2 -")/ e . Approximating the time interval by the average of the distribution Pp(r), 
t « (t) = 1/v, we obtain the trajectory x{t) in the form of a sum of mutually independent variables with the same 
distribution 




(15) 



x{t) 



v(t) 



(16) 



n 




(17) 



where the n'th interval contains the time t. The distribution of the variable 77 = ^ 



l + (a-2)/6» 



> Pvi 7 !), is given by 
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FIG. 3: (Colour online) Probability density distributions of the variable x for the case of the exponential covariance at the time 
t = 5. The curves correspond to a = 0.8, 1.2, 1.5, 1.8 and 2 (from right to left); Vq = 1. The red dashed line results from 
inversion of the characteristic function (|21|) for a, = 1.5 and 7 = 1/a. 




FIG. 4: (Colour online) Time evolution of the density distribution p x (x,t) for algebraic C(t) and v(t) = £(t). The curves 
correspond to the following times: 5, 10, 15, 20, 25 and 30 (from left to right). Inset: variance as a function of time; the 
straight line marks the dependence t ' 83 . 



P n {v) = Pi[€(v)]\d£/dr)\ ~ The variance of the distribution P v (rf) is finite if 

a.6 



a + 6-2 



> 2. (18) 



If the condition (|18[) is satisfied, the variable x has finite variance as well and its distribution converges with time 
(n — > 00) to the normal distribution, though that convergence may be very slow. The variance is always finite for 
the case 6 = 1, whereas it is always infinite if 6 > 2. Therefore, a dynamical system in which the velocity is given 
by KP with the algebraic covariance can be characterised by the finite position variance, although the Levy stable 
distribution is assumed. This conclusion is the main result of the paper. 

The numerical results for the case v(t) = £(t), presented in Fig. 4, show that the tails of the probability distributions 
are exponential. The variance rises with time as t 83 up to t = 50, which dependence indicates subdiffusion. The case 
of stationary tails, v(t) = g(t)t;(t), is presented in Fig. 5. The results are qualitatively similar to the preceding case but 
the distributions are broader, and they expand faster with time. The shape of the tails is exponential, ~ exp(— 7|jc|), 
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FIG. 5: (Colour online) The same as Fig. 4 but with v(t) = g(t)£(t). The straight line in the inset marks the dependence t 1,65 . 

and the parameter 7 rises with time. The time dependence of the variance is governed by the relation t 165 up to 
t = 40, whereas for larger times the growth becomes even faster. Anyway, the dependence is stronger than linear. 

The Langevin equation with uncorrelated noise may predict the finite variance if that noise is multiplicative. That 
problem was discussed in Ref.(45| for the algebraic multiplicative factor |v| _e / Q , where 6 + a > 0. The conclusions 
depend on a particular interpretation of the stochastic integral. In the Ito interpretation, for which the corresponding 
Fokker-Planck equation contains a variable diffusion coefficient, the tails of the distribution are the same as for the 
input noise: On the other hand, the Stratonovich version predicts the tails and the process for 

the case without potential is subdiffusive if a + 9 > 2. Heavy tails with finite variance may also emerge when one 
introduces a deterministic force, both linear and nonlinear (46. \i7\. 



IV. COMPARISON WITH THE ORNSTEIN-UHLENBECK-LEVY PROCESS 

The Ornstein-Uhlenbeck process describes Brownian particle velocity in the framework of the standard Brownian 
motion theory [48[ . The process is defined by the Langevin equation 

v(t) = -jv(t) + f)(t), (19) 

where the white noise r)(t) generates a process with stationary and independent increments. For the Gaussian case, the 
covariance assumes the exponential form with the decay rate 7 and (flT?)) is the only process with the above properties, 
according to the Doob theorem [4{|. In the case of the general stable distributions, Eq. (fl"9|) defines the Ornstein- 
Uhlcnbeck-Levy process. Its covariance is also exponential when understood in a sense of truncated distributions [50J. 
The probability distribution of the variable v(t) is given by the characteristic function [5i| 



p(k, t) = exp 



_J_|/t|«(l_e- Q ^) 

017 



(20) 



it converges with time to the stationary distribution. The two-dimensional process (x, v) is Markovian and the 
probability density p(x,v,t) determines the position distribution after performing integration over velocity [50j . Its 
characteristic function reads 

&(M) = e-'WI*l a , (21) 

where 



alt) = / dn (22) 

Jo 1 - « 

and g = 1 — e -7 *. According to Eq.(|20p. OUL process has the same time-scale as KP for v = 7a, both processes can 
then be compared. In order to invert the transform (|21l) we utilise the fact that p x {%, t) is the stable Levy distribution 
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FIG. 6: (Colour online) Density distributions p x (x,t) for v(t) which is constructed from the OUL processes according to 
Eq. (|23|) . The solid lines correspond to a — 1 , 1.2, 1 .5 and f .8 (from right to left). The solid straight lines mark the dependence 
a; -1-a! and the dashed green line the result for v(t) from KP with a — 1.5. All distributions were calculated at t = 5. Inset: 
autocorrelation functions calculated from the trajectory simulations with 1/1(7) = I/70 (lower curve) and 0(7) = 1/(2a/77o); 
in both cases 70 = 10. The lower and upper red dashed lines correspond to the functions (1 — e -7 ° r )/(7ot) and 0.15/^/r, 
respectively. 



and it can be expressed as the Fox function [52j | ; its numerical values follow from series expansions both for small and 
large x. A result of the expansions for a = 1.5 is compared with the corresponding distribution for KP in Fig. 3; both 
distributions agree. 

Processes with long tails of the autocovariance can also be constructed by OUL, namely as their superposition. 
Following Ref. [53j , we define a non-Markovian compound process in the form 

1 " 

"(*) = "!>(*)> (23) 

i=l 

where the atomic parts are OUL processes, 

«»(*) = -7i«i(*) + »)(*). (24) 

and n — > 00. The parameter 7 in the atomic process is a stochastic variable with a given probability distribution ^(7) 
which satisfies the normalisation condition J °° 0(7)^7 = 1. General solution of Eq. (j2"3")l and Eq. (pM|) over the entire 
real line is given by 

v(t)= lim / - y%-»<*-* >ri(fW. (25) 

Covariance follows from the formula for the atomic process, Eq.([5]), where the joint probability distribution has to be 
conditioned by ^(7). The final expression involves the Laplace transform: 

/>oo 

C(t) - / e-TtMdry. (26) 
Jo 

Eq. (l26[) can be inverted to obtain 4>(l) f° r a given form of the covariance. Two examples of the power-law form 
are presented in Fig. 6; C(t) was calculated by averaging over an ensemble of trajectories from Eq. (l24]) . The case 
C(t) = (1 — e _7ot )/(7ot) corresponds to 4>(l) — I/70 f° r 7 < 7o and zero elsewhere. The second example involves a 
slowly decaying tail, ~ i -1 / 2 , which appears when we choose ^(7) = 1/(2^/770) (7 < 70). The emergence of long tails 
in the autocorrelation function results from the presence of small values of 7 in the integral (|26p . 

Probability distributions of position, p x (x,t), for v(t) calculated from Eq. ([23|) . are presented in Fig. 6. All the 
curves, which correspond to different values of a, assume the power-law shape of the tails, la;! - " -1 . Therefore, the 
variance is divergent. Those results are qualitatively different from the distribution for the jumping process which is 
also presented in the figure. 
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V. SUMMARY AND CONCLUSIONS 

We have considered a Markovian jumping process, KP, which is characterised by long jumps and finite correlation 
time. Its key property is the dependence of the jumping rate on the process value which results, in particular, in a 
variety of correlation functions. The process value distribution is such that the corresponding master equation has a 
stationary solution in the form of a Levy stable distribution with divergent variance. The correlation functions can 
be properly defined if one introduces a cut-off at some large process value. Numerical analysis, performed for the 
stability index a = 3/2 and the power-law correlation function, shows that the stationary solution Pi(£) is unstable: 
it decays with time to the delta function. However, the shape of the tails converges with time to that of Pi(£) and 
a simple modification of the process - which consists in multiplying the process value by an appropriate function of 
time - makes the tails stable. 

When the particle velocity is assumed as the Levy distributed jumping process the transport properties of 
the system are qualitatively dependent on a particular form of its autocorrelation function. The tails of the position 
distribution for the exponential correlations have the shape the same as for the velocity, indicating infinite 

variance. It may not be the case for the power- law form of C(t) since long jumps are penalised by small time intervals. 
It has been demonstrated that if the parameters a and 9 satisfy the condition (fTS)) . the position distribution falls 
relatively fast making the variance finite. Suppressing of long tails of the distribution resembles the Levy walk 
approach where the jump length is governed by the time needed to perform the jump. The diffusion process for the 
case of stable tails appears anomalous; it is characterised by the variance which rises faster than linearly with time. 
On the other hand, diffusion may be accelerated if condition ([T8j) is not satisfied. 

The existence of finite variance of the position distributions - despite long velocity tails - is not a generic property 
of power-law correlated processes but rather a specific feature of KP. Other processes may exhibit different behaviour; 
we have demonstrated that a superposition of the OUL processes is characterised by tails of the form |:r| _a_1 and 
diffusion is accelerated, similarly to the case of the exponential correlations. Divergence of the variance does not 
violate physical princi ples for some problems, e.g. the diffusion on a polymer chain in chemical space or single 
molecule spectroscopy [17] . However, in dynamical problems with finite particle mass, a finite propagation velocity is 
required. From that point of view, KP offers a reasonable alternative to OUL as a model of the noise. 

APPENDIX 

Random numbers distributed according to are generated by inversion of the distribution function $(£). The 
following algorithm is restricted to the case a = 1.5 and 9 — 1. At first, we numerically evaluate z^(£) from Eq. (|10[) 
for consecutive values of £ up to !; m ax, with a small step A£, and store the results. ^ m ax means the process value 
for which the Levy asymptotics is already reached: -Pi(£) ~ £~ 2 - 5 > £ max ). Therefore, we have either v(£) = 

<■*,■■> n ax. ) or = v{^ m ax) + c(V? — V '6nai) (£ > Cmax), where c is a constant. Then we evaluate 

and Q(0 from Eq.@. We must invert the function $(£) = v{% )Pi(£')*:7(K£)) = r > where r € (0,0.5) is 
a uniformly distributed random number. $(£) is calculated numerically and stored if r < $(^ mol ). The asymptotics 
must be treated separately. Performing the integral yields the equation 

HUa X ) + cir 3/2 + car 1 + c 3 - r = 0, (Al) 

which resolves itself to a third-order algebraic equation. Finally, the sign is sampled. In the calculations presented in 
the paper, the constants have the following values: £ max = 20, c = 1.228, c\ — —2.162, c-i — —0.725 and C3 = 0.06038. 

The above algorithm may be generalised to other values of a and 9 but in the most cases Eq. (|A1[) must be solved 
numerically. Moreover, the integral in the expression for $(£) (£ > £ m ax) may not be an elementary function if 9 =/= 1. 
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